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Abstract. Recent observations of Cepheids using infrared interferometry and Spitzer photometry have detected the presence 
of circumstellar envelopes (CSE) of dust and it has been hypothesized that the CSE's are due to dust forming in a Cepheid 
wind. Here we use a modified Castor, Abbott & Klein formalism to produce a Cepheid wind, and this is used to estimate the 
contribution of mass loss to the Cepheid mass discrepancy Furthermore, we test the OGLE-III Classical Cepheids using the 
IR fluxes from the SAGE survey to determine if Large Magellanic Cloud Cepheids have CSE's. It is found that IR excess is a 
common phenomenon for LMC Cepheids and that the resulting mass-loss rates can explain at least a fraction of the Cepheid 
mass discrepancy, depending on the assumed dust-to-gas ratio in the wind. 
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INTRODUCTION 

Cepheid variable stars have been studied for more than 
200 years, since the discovery of variability in 8 Cephei 
|Q]] . Even now, there are still many unanswered questions 
about the structure and evolution of these stars and new 
questions are being raised. One such challenge is the dis- 
covery of circumstellar envelopes surrounding Galactic 
Cepheids using infrared interferometry |2|, |3l |4J . 

These envelopes appear as an IR excess which has 
been observed before. Galactic Cepheids have been ob- 
served using IRAS and IR excess was detected in some 
Cepheids flllla]. More recently, Spitzer observations have 
confirmed IR excess around Galactic Cepheids [7] . There 
is evidence for IR excess about I Car and RS Pup based 
on observations of the spectral energy distribution from 
V-band to 100 jJ.m [8, also see Alexandre Gallenne's 
poster]. These observations provide compelling evidence 
for the existence of circumstellar envelopes surrounding 
Cepheids. 

While observational evidence is mounting for the exis- 
tence of the circumstellar envelopes, there are a number 
of possible explanations for what they are. One possi- 
bility is that the envelopes are not envelopes at all but 
are instead disks J^]. This model was suggested as an 
alternative model for computing the distance to RS Pup 
from light echoes instead of the spherically symmetric 
model [ 10]. The disk model was proposed as an alterna- 
tive solution only and was not intended to necessarily be 
a physical explanation, however, it is important to con- 
sider the model. One argument against the disk model 



is the fact that interferometry has detected shells about 
every Cepheid that has been observed. If the IR excess 
were due to dusty disks then one would expect that the 
amount of excess would depend on the inclination of the 
disk which is random with respect to the line of sight of 
the observer. Hence one would expect only a small frac- 
tion of Cepheids to exhibit an IR excess, contrary to what 
has been observed. 

Another possibility is that the envelopes are relics 
from earlier stages of stellar evolution. For instance, it is 
argued that the cold component of the nebula surround- 
ing RS Pup is interstellar in nature [8], and the argument 
is supported by observations of a circumstellar envelope 
surrounding the Herbig Be star HD 200775. The Be star 
may be an analog to an earlier stage of evolution for 
RS Pup. The connection is strengthened by the observed 
elongation of the nebula, suggesting bipolar mass loss 
fiUl . However, there is no evidence that the warm cir- 
cumstellar envelopes are evolutionary relics, and inter- 
ferometric observations of the non-pulsating yellow su- 
pergiant a Persei suggest that the CSE's are related to 
the evolution of Cepheids. 

A simple explanation for the existence of CSE's is 
mass loss from the Cepheids themselves. This is not a 
new idea, Cepheid mass loss has been modeled and it was 
found that Cepheids could lose about 10~ 7 M@/yr. This 
amount of mass loss would have a significant effect on 
the evolution of a Cepheid [12]. The concept of Cepheid 
mass loss has been "rediscovered" with the observations 
of CSE's. Furthermore, spectroscopic observations of 
a number Cepheids found asymmetry in the Ha line 



consistent with an outflow iTHl . 

If the CSE's are generated from Cepheid mass loss 
then this mass loss affects the IR Leavitt Law (Period- 
Luminosity relation), may add an additional uncertainty 
to the IR surface technique for determining the angular 
diameter of Cepheids, and play a role in the evolution 
of these stars. It is important to characterize the nature 
of the CSE's and the mass-loss mechanism of Cepheids. 
One approach to understanding CSE's is to model the IR 
excess of a large sample of Cepheids. Fortunately, there 
exists IR observations of a large number of Cepheids 
thanks to the correlation of the OGLE-III Cepheids 13 
with the SAGE survey of the Large Magellanic Cloud 
iflHl as well as the 2MASS survey G3l- This correla- 
tion yields VIJHK and IRAC magnitudes for about 1800 
Cepheids for each of the two epochs of SAGE observa- 
tion plus the average flux from the two epochs. Using this 
data we can test the existence of CSE's in LMC Cepheids 
and if the source of the CSE's is mass loss. 



METHOD 

In this work, we test the existence of CSE's surrounding 
LMC Cepheids by fitting the spectral energy distribution 
of each Cepheid in the sample. The brightness as a func- 
tion of wavelength is fit assuming a stellar luminosity 
plus a dust shell luminosity, where the dust forms at a 
larger distance from a Cepheid and is due to mass loss. 
The stellar luminosity is 

L v (Star)=47CR 2 7CB v (T eS ), (1) 

and the shell luminosity is 

L v {Shell) = ^ <a \ > ^< 

4-71 < fl 3 > pv dm 

xQ^ T B v {T dmt )[\-W{r)]dr. (2) 

JR C 

The shell luminosity depends on the mean surface area 
and volume of the dust grains as well as the mean density 
of the dust grain. The dust is assumed to be silicate dust 
with a mean density of p = 3.7 g/crrr 1 , the dust ranges in 
size from a = 0.005 to 0.25 jJ.m and we use the Mathis et 
al. dust size distribution B17I1 for computing < a 2 > and 
< a 3 >. The dust velocity is about 100 km/s and the ab- 
sorption coefficient Q§ °< 1/A for A < 10 pirn. The term 
W(r) is the dilution factor and is proportional to r~ 1/2 . 
The inner boundary is the condensation radius, defined 
as the distance from the star where the gas reaches a 
temperature T — 1200 K. The dust is assumed to be in 
equilibrium with the stellar radiation and hence the dust 
temperature T dust = T eS W{r) 1 / 4 d. 

The stellar temperature is estimated using a 
temperature-color-period relation for LMC Cepheids 



lfl9tl . Thus if we can guess the radius and the mass-loss 
rate of a Cepheid then we can compute the total lumi- 
nosity (sum of stellar plus shell) and assuming a distance 
modulus of 18.5, we predict the apparent magnitudes 
of each Cepheid as a function of wavelength. Therefore 
we X 2 -fit these two parameters for each Cepheid in the 
sample and predict radii and mass-loss rates. 

RESULTS 

The values of the x 2 fit f° r eac h epoch of SAGE obser- 
vation is shown in Figure [TJ Because of the large spread 
of values we apply a ^ 2 -cut of 1.25 to the data. Also, one 
might argue that the because of the large spread of % 2 
values that mass loss is not a reasonable model but we 
have computed % 2 fits of radius only and find a similar 
spread of fits. The predicted radii and mass-loss rates of 
the Cepheids with % 2 < 1.25 are shown in FiguresQ]and 
|2] It is worth noting that fitting a Period-Radius (P-R) re- 
lation to the predictions yields a slope of 0.68 consistent 
with observed P-R relations. The mass-loss rates range 
from about 10~ n to 10~ 7 M & /yr suggesting that mass 
loss is an important phenomena in Cepheids. 

One question is what is the driving mechanism for 
mass loss in Cepheids. In earlier works I20l 12111 . the 
mechanism was explored using a sample of nearby 
Galactic Cepheids Il22ll . The Castor, Abbot and Klein 
(CAK) model for radiative-driven winds 02311 is applied 
to the sample and radiative-driven mass-loss rates are 
computed. These rates are shown in Figure[3]and a visual 
comparison of the predicted radiative-driven mass-loss 
rates with the best-fit rates for the LMC Cepheids sug- 
gest that the mass loss cannot be due to radiative driving. 
The mass-loss rates for short-period Cepheids are signif- 
icantly smaller than the computed rates and there is an 
obvious correlation between the radiative-driven mass- 
loss rates and the pulsation period that is not seen in the 
LMC Cepheids. 

It can easily be concluded that radiative driving is an 
insufficient driving mechanism. We have proposed an 
alternate theory that mass loss is driven by pulsation 
and shocks generated in the atmosphere due to pulsation 
similar to earlier works [12]. The idea is tested using 
modified version of the CAK method to include the effect 
of pulsation and shocks. The shock velocities are taken 
from the hydrodynamic model of 8 Cephei B2411 with 
a scaling relation derived to compute shock velocities 
for other Cepheids. The pulsation/shock-driven mass- 
loss rates for the sample of Galactic Cepheids is shown in 
Figure|3]as a function of pulsation period. It is found that 
the pulsation/shock-driven mass-loss rates do not have a 
period dependence. The mass-loss rates also range from 
10~ 10 to 10~ 7 M Q /yr, although only a few Cepheids 
have large mass-loss rates. 
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FIGURE 1. (Left) The values of x 2 for the fit of the radius and mass-loss rate for each Cepheid in the OGLE-III sample for each 
epoch of SAGE observation as a function of period. The horizontal line represents the cut-off of % 2 = 1.25. (Right) The best-fit 
radius of each Cepheid where the value of y} < 1.25. 



We also compare predicted mass-loss rates from for 
theoretical models of Galactic, LMC, and SMC 
Cepheids 112 ill to the computed mass-loss rates of the 
LMC Cepheids. The pulsation/shock-driven mass-loss 
rates are determined to be up to 2 x 10~ 7 M Q jyr which is 
comparable to the mass-loss rates determined here from 
observations. This is further evidence that the mass-loss 
mechanism for Cepheids is due to pulsation and shocks. 

The pulsation/shock theory for Cepheid mass loss 
roughly agrees with behavior of the best-fit mass-loss 
rates for the LMC Cepheids as a function of period, 
but seems to predict smaller mass-loss rates. One rea- 
son the mass-loss rates may be underestimated is that the 
modified-CAK method requires knowledge of the mass 
of a Cepheid and this was done assuming a Period-Mass- 
Radius relation derived from adiabatic pulsation mod- 
els. It is likely these masses are overestimated and if the 
masses are actually smaller then the mass-loss rates may 
increase significantly. Another source of error is that the 
mass-loss model is a quasi-static model and ignores time- 
dependent dynamics in the atmosphere of a Cepheid. 



CONCLUSIONS 

In this presentation, we have shown computed mass- 
loss rates for LMC Cepheids from the OGLE-III survey 
that have been correlated with the 2MASS and SAGE 
surveys to determine IR fluxes. The mass-loss rates of 
Cepheids are found to range from l(T 10 to 1(T 7 M Q /yr. 



These mass-loss rates contribute a shell flux that is about 
10 - 30% of the stellar flux at 8.0 \im. 

These rates are not strictly period-dependent as seen 
for prediction of radiative-driving. There is no obvious 
reason to expect that Cepheid mass loss should be a func- 
tion of pulsation period, especially since a period depen- 
dence is equivalent to a Reimer's relation for Cepheid 
mass loss. A Reimer's relation suggests that the potential 
energy of the wind is proportional to the stellar luminos- 
ity, however, shocks in the atmosphere are not necessar- 
ily proportional to the luminosity alone. The shocks are 
found to be non-linearly proportional to the luminosity, 
pulsation period and mass 12011 . which leads to a signif- 
icant deviation from a Reimer's relation. Therefore it is 
reasonable for the mass-loss rates to not be a strict func- 
tion of period. 

Mass loss is a reasonable explanation for the existence 
of the CSE's surrounding Cepheids and hence suggests 
that there is an IR excess for the majority of Cepheids. 
This IR excess affects the structure of the IR Leavitt 
Law |Hl|2l,|27l|2i] by making Cepheids appear brighter 
than they actually are. This result is not significant if the 
amount of IR excess is similar for all Cepheids indepen- 
dent of metallicity but if the IR excess depends on metal- 
licity then mass loss will cause an additional uncertainty 
to the IR Leavitt Law for computing distances. 

If the mass loss affects the IR Leavitt Law, then it is 
also reasonable to expect that mass loss will contribute to 
the IR surface brightness technique and potentially cause 
Cepheid angular diameters to be overestimated. It is im- 
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FIGURE 2. The best-fit mass-loss rates for the LMC Cepheids from the OGLE-III survey with % 2 < 1.25. for each epoch of 
SAGE observation. 



portant to characterize the how much angular diameters 
are overestimated if one wishes to compute distances to 
Cepheids using the Baade-Wesselink method. 

Cepheid mass loss may also be a solution to the 
Cepheid mass discrepancy II29[ 13011 . where there are two 
arguments for the current value of the mass discrep- 
ancy. Caputo et al. ll29ll argues the mass discrepancy is 
about 20% for small mass Cepheids (M w 4 M Q ) and de- 
creases as a function of mass. On the other hand, Keller 
lf30h argues the mass discrepancy is about 17% for all 
masses. If we consider the analysis in this work, short- 
period Cepheids (logP < 1), with an average mass-loss 
rate of 10~ 7 M Q /yr will lose about 1 M over a cross- 
ing of the instability strip which is consistent with both 
arguments for the mass discrepancy. However, for long- 
period Cepheids the total mass lost over a crossing of 
the instability strip will a smaller fraction of the stellar 
mass because the evolutionary timescale decreases with 
increasing mass. That result is consistent with the argu- 
ments of Caputo et al. but not of Keller. It is necessary to 
better understand both Cepheid mass loss and the mass 



discrepancy to determine if mass loss is truly the solution 
for the mass discrepancy. 

REFERENCES 

1. J. Goodricke, and Bayer, Royal Society of London 
Philosophical Transactions Series I 76, 48-61 (1786). 

2. P. Kervella, A. Merand, G. Perrin, and V. Coude Du 
Foresto, A&A 448, 623-631 (2006). 

3. A. Merand, P. Kervella, V. Coude Du Foresto, G. Perrin, 
S. T. Ridgway, J. P. Aufdenberg, T. A. Ten Brummelaar, 
H. A. McAlister, L. Sturmann, J. Sturmann, N. H. Turner, 
and D. H. Berger, A&A 453, 155-162 (2006). 

4. A. Merand, J. P. Aufdenberg, P. Kervella, V. C. d. Foresto, 
T. A. ten Brummelaar, H. A. McAlister, L. Sturmann, 

J. Sturmann, and N. H. Turner, ApJ 664, 1093-1101 
(2007). 

5. C. W. McAlary, and D. L. Welch, AJ 91, 1209-1220 
(1986). 

6. H. P. Deasy, MNRAS 231, 673-694 (1988). 

7. M. Marengo, N. R. Evans, P. Barmby, G. Bono, and 
D. Welch, "Spitzer Detection of Cepheid Circumstellar 
Emission," in The Evolving ISM in the Milky Way and 



10" B 



10- 7 ; 




1 10 100 1 10 100 

Period (d) Period (d) 



FIGURE 3. (Left) The predicted radiative-driven mass-loss rates for Galactic Cepheids computed using the CAK method. (Right) 
The predicted pulsation/shock-driven mass-loss rates for the same sample of Cepheids. The dotted line represents the relation 
between the radiative-driven mass-loss rates and the pulsation period. 



Nearby Galaxies, 2009. 

8. P. Kervella, A. Merand, and A. Gallenne, A&A 498, 
425-443 (2009). 

9. M. W. Feast, MNRAS 387, L33-L35 (2008). 

10. P. Kervella, A. Merand, L. Szabados, P. Fouque, 

D. Bersier, E. Pompei, and G. Perrin, A&A 480, 167-178 

(2008) . 

11. H. E. Bond, and W. B. Sparks, A&A 495, 371-377 (2009). 

12. W. M. Brunish, and L. A. Willson, "Evolution of 
Cepheids with Pulsationally Driven Mass Loss," in Stellar 
Pulsation, edited by A. N. Cox, W. M. Sparks, and S. G. 
Starrfield, 1987, vol. 274 of Lecture Notes in Physics, 
Berlin Springer Verlag, pp. 27-+. 

13. N. Nardetto, J. H. Groh, S. Kraus, F. Millour, and 
D. Gillet, A&A 489, 1263-1269 (2008). 

14. I. Soszynski, R. Poleski, A. Udalski, M. K. Szymanski, 
M. Kubiak, G. Pietrzynski, L. Wyrzykowski, 

O. Szewczyk, and K. Ulaczyk, Acta Astronomica 
58, 163-185 (2008). 

15. M. Meixner, K. D. Gordon, R. Indebetouw, J. L. Hora, 
B. Whitney, R. Blum, W. Reach, J.-P. Bernard, M. Meade, 
B. Babler, C. W. Engelbracht, B.-Q. For, K. Misselt, 

U. Vijh, C. Leitherer, M. Cohen, E. B. Churchwell, 
F. Boulanger, J. A. Frogel, Y. Fukui, J. Gallagher, 
V. Gorjian, J. Harris, D. Kelly, A. Kawamura, S. Kim, 
W. B. Latter, S. Madden, C. Markwick-Kemper, 
A. Mizuno, N. Mizuno, J. Mould, A. Nota, M. S. Oey, 
K. Olsen, T. Onishi, R. Paladini, N. Panagia, P. Perez- 
Gonzalez, H. Shibai, S. Sato, L. Smith, L. Staveley-Smith, 
A. G. G. M. Tielens, T. Ueta, S. van Dyk, K. Volk, 
M. Werner, and D. Zaritsky, AJ 132, 2268-2288 (2006). 

16. C.-C. Ngeow, S. M. Kanbur, H. R. Neilson, 

A. Nanthakumar, and J. Buonaccorsi, ApJ 693, 691-696 

(2009) . 

17. J. S. Mathis, W. Rumpl, and K. H. Nordsieck, ApJ 217, 
425-433 (1977). 

18. Z. Ivezic, and M. Elitzur, MNRAS 287, 799-8 1 1 (1997). 

19. J. P. Beaulieu, J. R. Buchler, and Z. Kollath, A&A 373, 
164-172 (2001). 

20. H. R. Neilson, and J. B. Lester, ApJ 684, 569-587 (2008). 



21. H. R. Neilson, and J. B. Lester, ApJ 690, 1829-1838 
(2009). 

22. P. Moskalik, and N. A. Gorynya, Acta Astronomica 55, 
247-260 (2005). 

23. J. I. Castor, D. C. Abbott, and R. I. Klein, ApJ 195, 
157-174 (1975). 

24. A. B. Fokin, D. Gillet, and M. G. Breitfellner, A&A 307, 
503-515 (1996). 

25. G. Bono, V. Castellani, and M. Marconi, ApJ 529, 
293-317 (2000). 

26. C. Ngeow, and S. M. Kanbur, ApJ 679, 76-85 (2008). 

27. W. L. Freedman, B. F. Madore, J. Rigby, S. E. Persson, 
and L. Sturch, ApJ 679, 71-75 (2008). 

28. B. F. Madore, W. L. Freedman, J. Rigby, S. E. Persson, 
L. Sturch, and V. Mager, ApJ 695, 988-995 (2009). 

29. F. Caputo, G. Bono, G. Fiorentino, M. Marconi, and 
I. Musella, ApJ 629, 1021-1033 (2005). 

30. S. C. Keller, ApJ 677, 483-487 (2008). 



